Two-Photon Transport in a Waveguide Coupled to a Cavity with a Two-level System 
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We consider a system where a waveguide is coupled to a cavity embedded with a two-level system 
(TLS), and study the effects when a two-photon quantum state is injected into the waveguide. The 
wave function of two outgoing photons is exactly solved using the Lehmann-Symanzik-Zimmermann 
(LSZ) reduction formalism. Our results explicitly exhibit the photon blockade effects in the strong 
atom-cavity coupling regime. The quantum statistical characters of the outgoing photons, including 
the photon bunching and anti-bunching behaviors, are also investigated in both the strong and weak 
coupling regimes. These results agree with the observations of recent experiments. 
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Introduction. — The Jaynes-Cummings (JC) system, 
involving cavity quantum electrodynamics (QED) for a 
two-level atom inside a cavity, is of most importance 
for quantum optics [l[ and its applications. In the past 
decade, the JC system, in the regime where there are only 
a few photons inside the cavity, has been very extensively 
studied due to its potential applications for both quan- 
tum information processing and quantum device physics 
• The later is usually based on some solid state systems 
resembling the JC system, such as the superconducting 
circuit QED systems Q and optomechanical architec- 
tures 0. Recent experiments about various JC systems 
include the observation of photon blockade in a macro- 
scopic cavity coupled to an atom Q , the demonstrations 
of on-chip cavities coupling to atom or atom-like objects 
as well as the integration of such atom-cavity sys- 
tems with waveguides [Zj. A few configurations, where 
the atom-cavity system is either side-coupled [3 HI, or 
directly coupled [j| [Io| to a waveguide, are schematically 
shown in Fig. QJ, and Fig. [T]d. 

In the context of these recent experiments, it is nec- 
essary to theoretically study the intensity and coherence 
properties of transmitted or reflected light, when a few- 
photon quantum state is injected into the system through 
a waveguide. In Ref. the quantum states inside 

the cavity were expanded on a basis of photon number 
states. Truncating the number of basis states then re- 
duced the Master equations to an ordinary differential 
equations which can be solved numerically. This system 
has also been simulated by using the quantum trajectory 
approach . Analytically, closed- form formulas regard- 
ing the transmission and coherence properties have been 
obtained, either in the weak excitation limit where the 
atom is assumed to be mostly in the ground state [l2| , or 
in a mean-field-like approach where the expectation value 
of operator product is taken as the product of operator 
expectation values [13|. 

In this Letter, using field-theoretic techniques 
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provide an exact analytic formula for the out-going pho- 
ton wave-functions, when a two-photon state is injected 



FIG. 1: (Color online) Two kinds of coupling structures 
and the schematic of energy spectrum. (a) The side- 
coupled waveguide; (b) The direct-coupled waveguide; (c) the 
schematic for the energy spectrum of JC model in the sub- 
spaces n = 0, 1, and 2. 



into the system shown in Fig. [T] Our work is in con- 
trast to most existing theoretical works that used coher- 
ent state input. Given that these systems will be ul- 
timately used to process quantum information, under- 
standing their response to states other than coherent 
states is very important. Also, our result is exact, with- 
out the need of either the mean-field approximation, or 
being restricted to the weak excitation limit. In this 
system, the indirect photon-photon interaction is strong 
when the atom has significant probability in the excited 
state jl3j . Moreover, it is known even a single-photon 
pulse, when properly designed, can completely invert an 
atom in this system (l5| . Thus, it is important to go 
beyond the weak excitation limit. 

Transport model and S -matrix. — The systems in Fig. 
QJi, with a waveguide side-coupled to a cavity, is described 
by a Hamiltonian [l6] H = Hyf + Hjq + Hi contain- 
ing three parts: (a) The waveguide Hamiltonian Hy? = 
Sfe £ fe a l a fe' where a& (at) denotes the annihilation (cre- 
ation) operators of the photon. = v \k\ is the waveg- 
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uide dispersion relation, and we take the speed of light v 
as unity, (b) The JC Hamiltonian for the coupling of the 
cavity field to the TLS 

Hjc = uj c a^a + fl |e) (e| + g{a) \g) (e| + a \e) (g\), (1) 

where a (a^) denotes the annihilation (creation) oper- 
ators of the photon in the cavity with frequency u c , 
|e) (\g)) denotes the excited (ground) state of the TLS 
with energy level spacing SI, and g is the coupling con- 
stant of the TLS and the cavity field, (c) The term 
H\ = V^2 k (a\a->rY{.c.)/\fL, describing the coupling be- 
tween the cavity and the waveguide. This term results in 
the decay of the cavity mode into the waveguide. Here, 
V is the waveguide-cavity coupling constant and L is the 
length of waveguide. 

In term of the bonding and anti-bonding waveguide 
mode operators defined as efc(ofc) = (a k ± a_j.)/v / 2, we 
can write H = H c + H Q . The bonding sub-system is 
described by 



where the bare Green function Go defined 
by iGq 1 (ui p ,p) = u p - s p + i0 + . 

For one or two-photons, this computation results in 
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where V = \[2V . The anti-bonding sub-system is de- 
scribed by H = J2k>o k°k°k and is decoupled from the 
cavity. 

Next, we utilize the Lehmann-Symanzik-Zimmermann 
reduction approach [14| to study the bonding sub-system 
as described by the Hamiltonian H c . We aim to calculate 
the S'-matrix elements between the incoming and outgo- 
ing n-photon states, specified by their photon momenta 
k = ki, k n and p = p±, ...,p n , respectively. The two- 
photon S'-matrix has the form 

Sp 1 p 2 k 1 k 2 — Sp 1 k 1 Sp 2 k 2 ~\~ Sp. 2 k 1 Sp 1 k 2 + iTp 1 p 2 k 1 k 2 7 (3) 

where S p k = 5 p k + iT pk is the single-photon S'-matrix 
element. To determine the T-matrix, we first exactly 
calculate the connected Green function G[ p; k] (w P , Wk) = 

/G [p;k] (t',t)nj 1 = i[-exp(ia;p j ^ - iu) kj tj)dt j dt' j /2'ir] in S-matrix element S* 



Here, G4 = CTafa^afa^a 'fa)a 'fa)) is the four point 
Green function of a(t) = exp(ii7 c fft)aexp(— iHcfft), and 
the effective non-Hcrmitian Hamiltonian H c g is obtained 
by simply replacing lu c with a = cu c — iV 2 /2 in -Hjc- 
The imaginary part of a is the cavity decay rate. {T ■■■) 
is the time-ordered average on the state |0) c \g). Here, 
|0) and \g) denote the vacuum state of cavity field and 
the ground state of TLS, respectively. Since excitation 
number N = a^a + |e) (e| commutes with H c g, in its in- 
variant subspace with N = n, H c h is diagonalized with 
the eigenstates |A„±) = N n ±{— y/ng \n — l) c |e) + [fi + 
(n — l)a — X n ±) \n) c \g)} and the corresponding eigenval- 
ues A„± = {Vt + (2n-l)a±[{VL- af+Ang 2 ] 1 / 2 } /2, where 
Af n ± are the normalization constants. The schematic for 
the spectrum of JC model is shown in Fig. [TJ;. we notice 
that the eigenstates \X n ±) are not orthogonal to each 
other. Thus, in the following calculations, we need to 
use the bi-orthogonal basis approach [17] with the eigen- 
states |A* ± ) of H* s corresponding to the eigenvalues 



The orthogonal relations are (A* T |A„±) = and 
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<A* ± |A. 

In the above bi-orthogonal basis, we can evaluate the 
correlations (l~afa)a)fa) and G4, to obtain the single- 
photon and two-photon T-matrices and S-matrices. The 
results are listed as follows. 1. For the bonding modes, 
the single photon S-matrix is S p k = tk5k P , where t k = 
exp(— i25k) and phase shift 5k — arg[(fc — Ai+)(k — Ai_)]. 
The anti-bonding modes are free of coupling, thus possess 
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the frequency domain with ui^ — u) kl , ...,u)k„ and w p = 



., 0J Pn , where 
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pk ~ ufe P' 

transmission coefficients are obtained as f k = fa — l)/2 
and tk = (tk + l)/2, which agree with Ref. [16|. 3. The 
two-photon T-matrix elements are explicitly obtained 



Vk=0 

n* k =o 
(4) 

with t = ti,...,t n , and t' = t[,...,t' n . Z[q k ^* k ] = 
J D[e,a,a]exp{i J dt[L e + ^kyoiVk^ + h.c)]} is a gen- 
erating functional with L c being the Lagrangian for the 
bonding subsystem, and a denoting the variables of the 
TLS. The T-matrix elements can then be related to the 
Green function by 
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where E = ki + k% is the total energy of the incident 
photons. The square norm T2 — |T PlP2 fe 1 fe 2 | 2 of T-matrix 
element exhibits the two-photon background fluorescence 
in the bonding mode. 

Two-photon wave functions — Let olr(xi) {cll,{x2)) 
denote the annihilation operators of right (left) moving 
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FIG. 2: (Color online) The two-photon background fluores- 
cence and the correlation functions for strong coupling regime 
g — 5, the system parameters are uj c — Q, — 10 and V is 
taken as unit: (a) Two-photon background fluorescence for 
E = (b) g (2) (0) of the reflected (solid (blue) curve) and 
transmitted (dashed (red) curve) photons; (c) g^(r) of the 
reflected (dashed (red) curve) photons for E/2 = Xf'J and 
transmitted (solid (blue) curve) photons for E/2 = SI; (d) 
(3 (2) (t)- 1 )/ 1 15 of the reflected (solid (blue) curve) photons 
for E/2 = fi and {g (2 \r) - 1)/10 6 of transmitted (dashed 
(red) curve) photons for E/2 = xfj. 



photons [IJ]. It follows from Eqs. © and © that the 
out-going state |X out ) = \t out ) + |r out ) + \rt oui ) for two 
incident right-moving photons with momenta fci and k 2 
contains three parts: (a) The quantum state of two trans- 
mitted photons 

I Wt) = / dxidx2t2(xi,x 2 )a\ { (xi)a' R (x2)\Q)\g) 



explicitly defined by the two-photon wavefunction 
t 2 (xi,x 2 ) = ^-e lEx -[t kl t k2 cos(A fc x) - F(X, x)], (8) 

Z7T 



where Afe = k\ — k 2 , 
F(X,x 



V 4 9 4 E s= ± S ( E - 2Ai.) exp[*(| - Ax,_ s ) \x[ 



explicitly defined by 
1 



(c) the left-right entangled two-photon state 

\rt out ) = / dx 1 dx 2 rt 2 {x 1 , x 2 )a j ' L (x 1 )a' f R (x2) |0) \g) 



describes the scenario where one photon is transmitted 
while the other is reflected, where 



rt 
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2^ 



e^ x [(t kl r k2 e 2lA ^+t k2 f kl e 



-2iAkX c 



2F(X,2x c )}. 
(11) 

Photon statistics can be studied through the coher- 
ence functions [1] 3 (2) (r) = G (2) (r)/ |G*W(0)| 2 , where 
G«(r) = (F ont \a F (x + r)a F (z) |F out ), G^(r) = 
(F ou t\ a F (x)a F (x + r)a F {x + r)a F (x) \F out ), and F = 
R and L correspond to the transmitted photons 
and reflected photons respectively, while |i? ut) = 
I tout) / (tout I tout) and | Lout) = |r ou t) / {rant kout) ■ We 



\t 2 (x,x 

|2 



for 



4(A 1+ - A x _) Yl s= ±[(E - A 2s ) n i= i, 2 (^ " Ai s )] ' 

(9) 

and x = x\ — x 2 and x c = (x± +x 2 )/2 are the relative and 
center of mass coordinates, respectively; (b) The quan- 
tum state of the two reflected photons 

I rout) = / dxidx 2 r 2 (xi,x 2 )a' L (xi)a' L (x2)\0)\g) 



r 2 {xi,x 2 ) 



iEx c 



2vr 



'f kl f k2 cos(A fe x) - F{X, x)]; (10) 



have 5 (2) (t) = C(t)/D, where C(r) 
photons in transmission, and \r 2 (x, x + t)|' s for photons 
in reflection, which is independent of x and D is the nor- 
malization constant. 

Strong coupling regime. - To explore the physical 
consequence of the result above, we first consider the 
strong coupling regime with g > V 2 . The two-photon 
background fluorescence T 2 and the correlation function 
g^(r) are shown in Fig. [2]on resonance, i.e., uj c = Q, and 
we assume the same energy for the two incident photons, 
i.e. A k = 0. When the average energy of two photon 
E/2 = X^±, T 2 has one sharp peak at A k — A p = 0, and 
the four peaks emerge (see Fig. [2^,) for E = X 2 °±- 

In Fig. [2Jd, g( 2 \0) as the functions of energy E/2 
per photon are shown for two reflected and transmitted 
photons, respectively. The two reflected photons exhibit 
sub-Poissonian statistics (g {2) {0) « 1) for E/2 = A^2, 
and super-Poissonian statistics (<7^(0) >> 1) for E/2 = 
X 2 °±/2 or Q. For E/2 = A^, the correlation functions 
g( 2 )(r) of two reflected (dashed (red) curve in Fig. O;) 
and transmitted (dashed (red) curve in Fig. f2Ji) photons 
exhibit anti-bunching and bunching behaviors, respec- 
tively. When E/2 = ft, g i2) (r) of two reflected (solid 
(blue) curve in Fig. [2ji) and transmitted (solid (blue) 
curve in Fig. [2];) photons exhibit large bunching and 
anti-bunching behaviors, respectively. Our results for the 
two reflected photons in the case of a side-coupled cav- 
ity exactly agree with that for the transmitted photons 
observed in the experiment [5[ for the case of a directly 
coupled cavity as shown in Fig. [TJa. This agreement 
is not accidental, since the transmitted photon states in 
the direct-coupled case can be mapped into the reflected 
photon states in the side-coupled case, using a transfor- 
mation detailed in Rcf. [lij . 

In the side-coupled structure considered here, the re- 
flected light arises purely from the decaying amplitudes 
from the cavity. Thus, single photon reflection peaks 
at a single photon energy of A^, which is the energy 
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FIG. 3: (Color online) The two-photon background fluores- 
cence and the correlation functions for weak coupling regime 
g — 0.5, the system parameters are the same as that in Fig. 
2: (a)-(c) represent the same physical meaning in Fig. 2; (d) 
(3 (2) M-l)/10 n of the reflected (solid (blue) curve) photons 
for E/2 = CI and (g {2) (r) - 1)/10 10 of transmitted (dashed 
(red) curve) photons for E/2 = XfJ. 



behaviors exhibited by g^ 2 \r) of two reflected photons 
for E/2 = A^°2 (dashed (red) curve in Fig. [3fc) and CI 
(solid (blue) curve in Fig. [5Ji) also illustrate the vanish- 
ing of photon blockade effect. In addition, g^ 2 \r) of anti- 
bunched and bunched transmitted photons for E/2 = CI 
and A ^2 are shown by the solid (blue) curve in Fig. [5J; 
and the dashed (red) curve in Fig. 

Conclusion. — We have analytically studied the two- 
photon transport in a waveguide coupled to the cavity 
containing TLS, and obtained an exact analytic formula 
describing the photon blockade effect. The quantum 
statistics of the outgoing photons are discussed in de- 
tails using the exact two-photon wave-functions. These 
results agree with the observations of recent experiment 
Our theoretical approach can also be generalized to 
the many-photon transport cases. 

The work is supported by National Natural Science 
Foundation of China a under Grant No. 10874091 . 



level for one-photon dressed state in the cavity. Simi- 
larly, two-photon reflection peaks when the two photon 
energy is at Ag±, where the cavity supports two-photon 
dressed states. However, since A^°2 ^ 2A^°2, two pho- 
tons each with energy E/2 = XfJ. is off resonance from 
the two-photon dressed state. In such a case, the sin- 
gle excitation by the first photon in fact prevents the 
second photon from entering the cavity, resulting in the 
photon-blockade effect. Therefore, in contrast to the case 
with direct coupling where the photon-blockade effect 
manifests as a vanishing two-photon transmission, in the 
side-cavity case the photon blockade effect manifests as a 
vanishing two-photon reflection effect. When EJ2 = CI, 
the single-photon reflection coefficient vanishes [16J . The 
two-photon reflection is due purely to the correlation in- 
duced by the TLS, which creates a two-photon bound 
state, and hence generates a large bunching effect [Icj . 
We, therefore, for the first time, provide an exact ana- 
lytic formula for the photon correlation function for this 
system, the special case of which agrees with the existing 
experimental data. 

Weak coupling regime. - For the weak coupling 
regime with g < V 2 , the two-photon background fluores- 
cence T2 and the correlation function g^ 2 '(r) are shown 
in Fig. [3J Fig. [3^, shows that the four peaks in Fig. [2jt 
merge into a single peak when E = X^l . 

We find that g& (0) of two reflected photons (Fig. [5b) 
has a simple structure with one peak at E/2 — CI and al- 
ways satisfies g^(0) > 1, which means that the statistics 
of reflected photons can not be sub-Poissonian and the 
photon blockade effect vanishes due to the small energy 
splitting of A^2 for the weak coupling g. The bunching 
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